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Abstract 

We here study the collective excitations of a dilute spin-polarized Bose-Fermi 
mixture at zero temperature, considering in particular the features arising 
from the interaction between the two species. We show that a propagating 
zero-sound mode is possible for the fermions even when they do not interact 
among themselves. 

PACS numbers: 03.75.Fi, 67.60.-g, 67.55.Jd 
Recent experimental progress in atomically trapped gases has led to a resurgence of 
interest in quantum fluids. A particular notable feature is the number of systems available, 
ranging from single component bose gas in the original experiments where Bose-Einstein 
Condensation (BEC) was first achieved |L| to binary bose mixture 0, spinor condensate in 
optical traps and degenerate fermi gas 0|| . Other systems have also received much recent 
attention, in particular bose-fermi mixture. This last mentioned system occurs naturally if 
'sympathetic cooling' is employed to reduce the kinetic energy of the fermions ||. There 
have already been several studies on the properties of this system. Questions addressed 
include stability against phase separation J7||| and collective excitations 0. 

Although Bose-Fermi mixtures have been studied intensively in low temperature physics 
in the context of 3 He- 4 He mixtures || , the atomically trapped gases offer many additional 
possibilities. By the choice of atoms, concentration of the various components, or the control 
of interaction strength among them by external fields ||10|| , one can unmask phenomena 



previously unobservable. In this paper, we shall study one example of this by considering 
the density oscillations of a bose-fermi mixture at low temperatures. We shall show that 
a variety of novel phenomena can arise due to the coupling between the two components 
for suitable parameters such as the ratio of the sound velocity of the bose gas to the fermi 
velocity of the fermions. In particular, we shall show that it is possible to have a propagating 
fermionic sound mode even in the absence of interaction among the fermions themselves. 
Sound propagation has also been considered in ref |J which however did not investigate the 
effects being studied here. We shall comment on this later. 

We shall then consider a mixture of weakly interacting bose and fermi gases at zero 
temperature. Both gases are assumed to be spin-polarized such as would be the case usually 
in magnetic traps. For a dilute mixture, interaction among the bosons themselves and 
between the bosons and fermions can be characterised by the scattering lengths and 
cibf in the s-wave channels. The fermions however do not interact among themselves since 
they are spin-polarized. For simplicity we shall consider a uniform system. We shall further 
assume that the gas is stable against phase-separation unless explicitly specified. We are 
interested in the density waves of this system. As we shall see in general the modes may be 
damped. Also since the density oscillations are likely to be studied by exciting the systems 
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with external potentials, we shall instead consider the density responses of the system under 
external perturbing potentials. Collective modes of the system will show up as resonances 
of these responses. 

The Hamiltonian density is given by 
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where the subscripts b and / denote bosons and fermions respectively, if)f, iftb are the field 
operators, m&, m/ the masses, fib, fJ-f are the chemical potentials, V b ext and Vf xt the external 
potentials. All ?/>'s and V ext 's are implicitly at the same physical point r in space. The 
interaction parameters g bb and ^ are related to the scattering lengths a bb and a b f by ^ = 
<iiih 2 a bb /m bb and = 2i\Ti 2 a b f / m r where m r is the reduced mass [ 
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We shall treat the interaction g bb and g b f within the Bogoliubov and random phase 
approximation respectively [JTT| . The results can be written in the physically transparent 
form: 



5n b (q,u) 
Sn f (q,u) 



-Xb[g bf Sn f + vn 

-Xf[gb f 5n b + V f ext } 



(2) 



expressing the response of the bosons and fermions to the potentials due to the other species 
and the external perturbations (the terms in the square brackets). Here 5n b (q,u)), 8nf(q,u) 
are the deviations of the bosonic and fermionic densities from equilibrium at wavevector q 
and frequency u, 
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and 
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are the (q and uj dependent) responses of the pure bosons and fermions systems respectively 
to effective external potentials. Nf = pfinf/2'K 2 is the density of states for the fermions. 
{jPf = (G^rif) 1 / 3 is the fermi momentum, Vf = pf/rrif) For simplicity, in eq (£|) I have 
already left out terms that are small if q « pf. uj should be interpreted as having a small 
and positive part. 

ec L@ can b e re-arranged as 
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Then finally 
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In the case where ^ = 0, 5rib = — XbV b cxt and 5nj = —XfVf Xt and the responses thus 
reduce to those of the pure bose and fermi gases. The corresponding formulas for \b and Xf 
were already given in eq (Bt) and eq (f§) above. Before we proceed we shall recall the behavior 



of these responses [1 T | and thus the collective modes. For simplicity we shall restrict ourselves 
to small wavevectors, i.e. q << m b Cb and pf, and without loss of generality uj > 0. The 
bosonic response Imx;, consists of a delta function at the excitation frequency uj = c b q. This 
is due to the Bogoliubov mode which is purely propagating and undamped. For the fermions 
however, there is no collective behavior. The absorptive part, Imxy, is finite for a whole 
range of frequencies |u;| < Vjq, known as the particle-hole continuum, arising from the many 
possibilities of independent particle-hole excitations. Re Xb is simple. It is given by g^ 1 at 
uj = and diverges to ±oo as uj —>■ c b q from below and above respectively. Re%/ is given 
by Nf at u — 0. It decreases with increasing uj, changes sign at around uj ~ 0.83v/q and 
approaches — oo as uj — > vjq from both above and below. For uj > vjq, it remains negative 
with its magnitude gradually approaching zero as u — » oo. 

Now we return to the bose-fermi mixture. The response 5rib to an external potential V b ext 
acting on the bosons only is given by Xbl (X~9bfXbXf)- The existence and the dispersion of the 
bosonic collective mode are determined by the solution to the equation (x?>) _1 — QbfXf = 0> 
i.e. 

[^+^+<£>v (£)(&■)*=<> (T) 

It will be convenient to discuss the normalized response 

Xb = 9bbXb/ (1 - glfXbXf) (8) 
Xb — 1 in the static limit (uj = 0, q —>■ 0) when there is no boson-fermion interaction 

Gfc/ = o). 

Similiarly the fermionic response to an external potential acting on the fermions alone 
is — QbfXbXf)- We shall discuss the behavior of the normalized quantity 

Xf = Nj l x f /(l-g 2 bf XbXf) (9) 

The normalization is chosen such that Xf = 1 m the static limit (u — 0, q — *■ 0) when there 
is no boson-fermion interaction (gbf = 0). 

Before proceeding let us first examine the responses at uj = 0. Stability requires that 
the density responses ~ 9bfXbXf) and — QbfXbXf) be positive. Using the uj = 

values of Xb and Xf above, these necessary conditions can be rewritten as > and 
W = Nfglf/gbb < 1- Using the expression of Nf given earlier, the last inequality gives 
n f9bf < \Agbb where A = ■^^.(^ 2 ) 2 ^ i as defined in 0. These conditions were derived 
earlier in J7| and || using slightly different considerations. For bosons and fermions with 
similar masses, we shall see shortly that W, a dimensionless parameter, serves as a useful 
measure of the coupling between the bosons and fermions. If the bosons and fermions have 
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similar masses, \W\ is of order \alf/abbTif 3 \ and thus typically small for dilute gases unless 
\a,bf \ >> \ a bb\- We shall limit ourselves only to the cases where |W|'s are small. 

We shall discuss now the behavior of Xb an d Xf i n turn. The results are qualitatively 
different depending on whether c^Vf. The velocity ratio u = c^/vf can be re-expressed 

as u = — ^^/e (jiba fl 1 ■ The value of u can basically be arbitrary without violating any 
stability criterion (not only the linear stability condition above but also others derived in 



Bosonic Response: 

1. Cb > Vf\ In this regime a propagating bosonic mode exists. It can be easily verified 
(e.g. graphically) that the mode frequency uj satisfies uj > Cbq ( > vjq). The original 
bosonic mode at uj = Cbq is pushed upwards by the particle-hole 'modes' lying below. Some 
examples are shown in Fig |l|. This mode 'repulsion' is generally expected (c.f. coupled 
harmonic oscillators). It is however of interest to examine the microscopic nature of the 
mode. At the mode frequency both Xb and Xf are negative. Thus, e.g., if gbf > 0, 5n b and 
5rif are of the same sign (see eqn @). The repulsion between the two species provides the 
enhanced restoring force and oscillation frequency. This frequency shift is typically small 
since usually W << 1. 

2. Cb < vf. In this case the original bosonic mode lies inside the particle-hole contin- 
uum of the fermions. The bosonic mode is thus Landau damped. For weak-coupling 
the damping, thus the width of the response, can be estimated easily using eq(^) to be 

~* N f9 2 bf ~ 



j (q><?)- Examples are shown in Fig ^. There is a small shift of the mode 

due to Rex/- The shift is towards higher frequency for u sufficiently close to 1 but opposite 
otherwise (Re^f < (>) for uj/vfq > (<) 0.83. ) 

3. It is also of interest to study the bosonic mode for g^ < 0. This is in fact the case for the 
6 Li- 7 Li mixture investigated in ref || , where the 7 Li bosons have a negative scattering length 
of ~ — 1.5nm. In this case the original bosonic system is unstable, and the Bogoliubov mode 
has an imaginary frequency for sufficiently small wavevector ( q < q c = 2mb\cb\/h, here 
\cb\ = [\gbb\ n b/ m b] 1 ^ 2 )- Since Njg^j > > g^, the system is still unstable in the presence of 
the fermions (see also above). Of interest is the effect of the fermions on the unstable 
mode. Now for imaginary frequencies u = ia, Xf(q^ a ) — ~ ^ [f ~ ^ an ~ 1 ^]} ^ s 

purely real and positive. Xf decreases monotonically with a from Xf — Nf at a — to as 
a — > oo. It can be easily verified that there is a real solution for a to the dispersion relation 
(cf. eq © ) 
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for sufficiently small q (which includes the physically most relevant region where a attains 
its maximum, i.e., the fastest growing instability). Thus the instability is not damped by 
the particle-hole degree of freedom. In fact it can be verified easily that, for given q, a is 
increased in the presence of the fermions. The system has become even more unstable. This 
mode has 5n h and 5rif of opposite signs and corresponds to phase-separation as expected. 
Fermionic Response: 

1. Cb > Vf\ In this case the fermionic response for < uj < qvf is only slightly modified. 
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A new feature appears near u ~ c^q > Vjq due to the coupling to the bosonic mode. An 
example is as shown in Fig |3[ 

2. c& < Vf\ In this regime there are two important features of the fermionic response. If 
u = is sufficiently close to 1, the imaginary part contains a sharp resonance at to above 
the particle- hole continuum (Fig|j). There are two ways of understanding this mode. It can 
be regarded as the continuation of the situation from c& > Vf, i.e., it is due to the bosonic 
mode which is itself slightly pushed up in frequency (c./., Fig |3|, note in particular the result 
for it = 1). Alternatively, this mode can be considered as a zero-sound mode induced by 
the bosons. The form for \f m ec L ® * s precisely that of an interacting fermi gas with 
s-wave interaction gjf (and therefore necessary with more than one spin species, where the 
response is given by Xf / (^+9ffXf)) though with an effective frequency dependent interaction 

9ff ~ *■ ~9bfXb, an effective s-wave Landau parameter given by F — > — 1]- The 

bosonic mode uj ~ c^q for q, sufficiently close to but below Vf will thus induce a zero-sound 



mode for the fermions just like an interaction among the fermions will. [12] Note however 
there cannot be a real s-wave interaction among the fermions as they are of equal spins. Thus 
this mode cannot be obtained by considering the effective interaction among the fermions 
as in ref 0. 

The frequency of this propagating mode can be estimated by using the well-known disper- 
sion relation of the zeroth sound uj/vjq ~ l + 2e~ 2( - 1+1 ^' ) with the effective F — » W/[(-) 2 — 1] 
as suggested above. In order for the velocity of the mode to be say 1% above Vf, then q, has 
to be within around 7% of Vf if W = 0.1. This estimate agrees with the numerical results 
of Fig |. 

The second interesting feature is that near the original bosonic mode frequency uj ~ q,g, 
there is a reduction in the absorptive part Imxf (see Fig H). In fact Inr^/ — > as uj — > c^q. 
This, as well as the corresponding behavior of Re%/, can be seen easily mathematically from 
eq(^|) due to the resonance nature of Xb at this frequency. Physically this can be regarded as 
due to mode-mode repulsion - the bosonic mode has pushed away the particle-hole 'modes' 
near u ~ q,g. This feature is present even for small coupling W . A larger W mainly increases 
the width of this 'transparent' region. Thus in fact the frequency dependence of Re%/ is 
actually stronger for smaller W's. 

The energy absorption by the bose-fermi mixture from an external perturbation acting on 
the fermions is thus substantially reduced for frequencies within this 'transparent region'. 
The width of this region can be estimated by using the observation that the fermionic 
response is roughly reduced by the factor 1 + f° r these frequencies. For the 

fermionic response to be reduced to say 1/2 of its bare value, then \u — Cbq\/cbq < W/2. 
This estimate agrees very roughly with the numerical results in Fig. ^. 

In conclusion, I have investigated the collective modes of a Bose- Fermi mixture, and have 
shown that there is important mode- mode coupling effects, especially if Vf ~ c&. 

I thank David Edwards for helpful correspondences. This research was supported by the 
National Science Council of Taiwan under grant number 89-2112-M-001-105. 
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FIGURES 




FIG. 1. Dimensionless bosonic re- 
sponses Rexfe and ImXb for u = Cb/vf > 1, 
W = 0.01. Lines for the imaginary parts are 
decorated with circles. 
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FIG. 2. Same as Fig. 1 but for 
u = q,/vf < 1, W = 0.01. 
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FIG. 3. Dimensionless fermionic re- 
sponses Rex/ and Imx/ for u = Cb/vf > 1, 
W = 0.1. The imaginary parts (lines deco- 
rated with circles) contain the particle-hole 
continua u < vjq and sharp spikes at the 
bosonic mode frequencies. 
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FIG. 4. Dimensionless fermionic re- 
sponse Imx/ for u = Cb/vf < 1 showing 
the zero-sound modes induced by the bosons. 
W = 0.1. Also shown is Imxj for the pure 
fermi gas (<#,/ = 0) for comparison. 
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FIG. 5. Same as Fig || except that now 
Rex/ is shown. 
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FIG. 6. Imaginary part of dimensionless 
fermionic response Im^/ for u = Cb/vf < 1, 
showing mainly the region u < vjq. 
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